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Abstract. We study the non-equilibrium relaxation of an elastic line described by 
the Edwards- Wilkinson equation. Although this model is the simplest representation 
of interface dynamics, we highlight that many (not though all) important aspects of the 
non-equilibrium relaxation of elastic manifolds are already present in such quadratic 
and clean systems. We analyze in detail the aging behaviour of several two-times 
averaged and fluctuating observables taking into account finite-size effects and the 
crossover to the stationary and equilibrium regimes. We start by investigating the 
structure factor and extracting from its decay a growing correlation length. We present 
the full two-times and size dependence of the interface roughness and we generalize the 
Family- Vicsek scaling form to non-equilibrium situations. We compute the incoherent 
scattering function and we compare it to the one measured in other glassy systems. 
We analyse the response functions, the violation of the fluctuation-dissipation theorem 
in the aging regime, and its crossover to the equilibrium relation in the stationary 
regime. Finally, we study the out-of-equilibrium fluctuations of the previously studied 
two-times functions and we characterize the scaling properties of their probability 
distribution functions. Our results allow us to obtain new insights into other glassy 
problems such as the aging behavior in colloidal glasses and vortex glasses. 

Keywords: Slow dynamics and aging (theory), self-affine roughness (theory), kinetic 
growth processes (theory) 
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1. Introduction 
1.1. Aim 

Interfaces are important in many physical, chemical and biological phenomena. In 
the physical context they appear in stochastic surface growth domain growth and 
coarsening phenomena (as domain walls) type-II superconductivity (as magnetic 
flux lines) [3], fracture cracks [4j, and fluid invasion in porous media |5j, among others 
realizations. Although specific problems involve different levels of complexity, as the 
inclusion of non-linear or disorder contributions, the main features of interface dynamics 
are already contained in the simplest theoretical formulation of the problem: the 
Edwards- Wilkinson (EW) equation |6]. 

The aim of this article is to exhibit, in a concrete and fully solvable example, 
several generic properties of the averaged and fluctuating aging dynamics of finite and 
infinite elastic manifolds. The problem we analyze is the relaxation of an elastic line 
described by the Edwards- Wilkinson equation [6]. The embedding space dimension 
plays an important role. In the particular one dimensional transverse space case on 
which we focus here the dynamics has aspects of diffusion, glassiness and saturation 
depending on the time and length scales observed. 

The relaxation process we are interested in is the following [7]. After equilibration 
at a temperature Tq the system is suddenly taken to a different working temperature T 
that can be higher or lower than Tq. Time is then set to zero. The line subsequently 
tries to adapt to the working temperature. The relaxation is characterized by the time- 
dependence of correlation and linear response functions. One lets the line relax until a 
waiting-time, t^, when the quantities of interest are first recorded and later compared 
to their values at a subsequent time t. 

A series of glassy properties of elastic lines evolving in disordered environments 
have been recently reported in different contexts, related to directed polymers in random 
media [SI [9], the vortex glass dynamics in high temperature superconductors p^[Tn[T2] . 
and domain wall motion in magnetic systems [131 E]- All these studies focused 
on models including disorder, which non-disordered counterparts belong to the EW 
universality class of interface dynamics [1], except for Ref. [14] in which the analysis 
is focused on the non-linear contribution in the Kardar-Parisi-Zhang (KPZ) [T5] 
universality class. In order to disentangle the effects of elasticity and disorder it is 
then important to study in detail the dynamics of the non-disordered counterparts. 

In general, the glassy phenomenon in finite dimensional models of elastic lines with 
and without quenched disorder appears as a dynamic crossover P [9l . For all waiting- 
times, t^, that are longer than a size, L, and eventually also temperature, T, dependent 
equilibration time, ti, the dynamics is stationary. Instead, for t^ < ti the system is in 
the preasymptotic regime and the relaxation occurs out of equilibrium as demonstrated 
by two-times correlations and linear responses that age. For each waiting-time the 
dependence on the time-delay. At = t — t^, then shows a growth regime for At < ti 
and saturation regime at longer time. At > t^, where the correlation functions saturate 



Out- of- equilibrium Edwards- Wilkinson equation 



3 



t ,+At 




z 



Figure 1. Two line configurations at different times, e.g. a:(z,tm) and x{z,tw + At), 
sliowing the relevant displacements defining the mean-squared-displacement of the line 
segments, B, the roughness, w^, and the mean-squared-displacement of the center of 
mass, D. 

to a size dependent value. 

Here we compute, analytically, the averaged two-times roughness and we relate 
it to the displacement field, which was the main focus of previous studies of the 
dynamics of elastic manifolds in random environments [U [TOl [16], and the center of 
mass evolution of the interface (see the sketch in Fig. [1]). The quadratic character of 
the model allows us to control the crossover to equilibration and saturation of finite 
lines. We next deduce several other correlators that are of interest in glassy dynamics. 
In particular, we analyze the incoherent scattering function and compare it to light 
scattering measurements in clay colloidal suspensions [17]. We interpret our results 
in terms of a two-times correlation length that we evaluate and confront to the one 
measured in other aging glassy systems [TH [T9|. We also study linear responses and 
their relation to the companion correlations. We discuss in detail the special features of 
these two-times observables linked to the multiplicative - diffusive - scaling form. 

Finally, we focus on thermally induced fluctuations. Recently, the importance of 
studying fluctuations - and not only averaged quantities - in dynamic phenomena was 
stressed in several contexts. Racz proposed to use scaling functions characterizing the 
fluctuations of global observables in non- equilibrium steady states to classify systems in 
'universality classes' dictated by symmetries and dynamic mechanisms [20]. In aging 
glassy systems the study of fluctuations seems to be fundamental to understand the 
mechanism for the dramatic slowing down and non-equilibrium relaxation. Chamon 
et al [211 [22] proposed a symmetry based sigma-model like theory for fluctuations in 
conventional glassy systems. For a number of reasons this theory is not expected 
to apply, without modiflcation, to interface dynamics. In more technical terms, the 
averaged interface dynamics is characterized by a multiplicative aging scaling that should 
result in the need to modify the approach in [21] to take this feature into account. 
We thus wish to confront the fluctuations of conventional glassy systems to those of 
interface models searching for similarities and differences. With this purpose we derive 
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the probability density functions (pdfs) of several two-times correlation functions. We 
compare our results to previous studies [23J, which focused on the time-delay dependence 
of the fluctuations 0, and to predictions of the time-reparametrization invariance scenario 
of glassy systems [HI [22] . 

The problem addressed here is related to a number of other models that have 
already been studied in the literature; among the pure cases one has the Langevin 
dynamics of the Gaussian and massless scalar field [23], the XY ferromagnet [211 [25] , 
the p = 2 spherical spin-glass [26] and the 0{N) ferromagnet in the large N limit |27j . 
Related models with quenched disorder are problems of elastic manifolds in quenched 
random environments [El [9l |13l HE]. In all these studies the finite size dependence was 
not taken into account and the dynamic fluctuations were not studied. We explain in 
the conclusions how our results relate to the ones in these papers. 



1.2. The model 

The EW equation for a scalar field x representing the height of a surface over a one 
dimensional substrate parametrized by the coordinate z (a one dimensional directed 
interface) is 

dtx{z, t) = vdlx{z, t) + f{z, t) + ^z, t), (1) 

2T 

t)) = 0, iaz, t)az', o) = —s{z - z')6{t - f), (2) 

7 

with u = c/j, c the elastic constant, 7 the friction coefficient, T the temperature of the 
thermal bath and (■ ■ ■) the average over the white noise ^, i.e. the thermal average. 
The term f = h/'y represents the effect of a perturbing field that couples linearly and 
locally to the height, —h{z, t)x{z, t). One can also consider other types of perturbations 
that couple to more complicated functions of the height, as discussed in Sect. [71 

We are interested in characterizing the dynamics of elastic lines with finite and 
infinite length L. Thus for convenience we shall take periodic boundary conditions in 
the ^-direction. Following Antal and Racz we introduce a Fourier representation of the 
position dependent distance of the line from its average 



6x{z, t) = x{z, t) - x{t) = J2 ^nit) e*^"^ (3) 

n=— 00 

with kn = 2TTn/L. The overline indicates an average over the full line's configuration; 

1 '■^ 



x{t) = — dz x{z,t) (4) 
Jo 

being then the mean height - or center of mass - of the interface. The Fourier coefficients 
are given by 

c„(t) = i [ dz [x{z,t) -xit)] e-'^"' (5) 
^ Jo 

I We note that the sohition presented in ^Sj corresponds actuaUy to the — solution in this article. 
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and their evolution is given by 

dtCnit) = -vklc^it) + Uit) + (6) 

2T 

{Ut)) = 0, {UtKn'it')) = —6^,-n'S{t - t') (7) 

7L 

for all n 7^ while Co(t) = at all t. Note that equation ([6]) effectively describes the 
position of a particle in a harmonic potential with spring constant vk"^ which implies an 
elastic constant softening with decreasing n or increasing system size L. The solution 
to equation ^ is 

c„(t) = c„(0) e-^^"* + e-^^"* /* dt' e'^^"*' [^^(tO + • (8) 

where we set the initial time to t = 0. The coefficients Cn(0) encode the structure of 
the initial conditions. We are interested in the evolution after an instantaneous quench 
from equilibrium at a generic temperature Tq to the working temperature T. We are 
then considering that the initial conditions satisfy 

(c„(0)c„(0)),,= |c„(0)p5„, (9) 

where |cn(0)p should reflect the equilibrium at Tq (see below). When the initial 
temperature is higher than the working temperature, Tq > T, the initial state is more 
disordered than the equilibrium one at T and the line tends to 'order' as time elapses. 
On the other hand, when Tq < T, temperature fluctuations roughen the initial - more 
ordered - configuration of the line. A special case is Tq = 0, corresponding to a perfectly 
fiat initial configuration. 

The EW energy is just the one of a massless scalar field and the model does not 
have a finite temperature static phase transition. 

Through out the article we shall present some figures which highlight the main 
analytical results. Without loss of generality we set 7 = z/ = 1 in all the figures. 
We distinguish noise averaged from fluctuating quantities by enclosing the former with 
angular brackets. In our expressions we do not write explicitly the T, Tq and L 
dependencies but one has to keep in mind that they are, in principle, always present. 

1.3. Organization of the paper 

The organization of the paper is the following. In Sect. 11.21 we had introduced the 
model and our conventions. Section [2] deals with the two-times structure factor and 
Sect. |3]with the associated two-times correlation length. In Sect. IHwe analyze the line's 
roughness and in Sect. [5] we study the displacement field as well as the motion of the 
center of mass. In Sect. [6] we derive a wave- vector dependent correlation inspired in 
the incoherent scattering function typical of particles in interaction. Section [7] analyzes 
different linear response functions and their relation to the associated correlations that 
is to say the modifications of the fluctuation-dissipation theorem (FDT). In Sect. [8] we 
study thermally induced fluctuations by computing the probability distribution function 
of the roughness, displacement and hnear responses and we relate to previous studies of 
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elastic lines in equilibrium as well as with the proposal for fluctuations in aging systems 
based on the development of time-reparametrization invariance. Finally, in Sect. [9] we 
summarize our results and we present our conclusions. The main features of the chosen 
observables are shown in flgures that display the analytic solution. 



2. Two-times structure factor 



A key quantity in the two-times evolution of the elastic line is the local displacement 
of the surface height, deflned as u{z,t,tw) = x{z,t) — which allows for the 

two-times generalization of several quantities. The structure factor associated to this 
displacement is 

{Sn){t,t^)=L{\Cr,{t)-cM\^) 



+ 



T 



(10) 

The structure factor deflned in this way is the two-times generalization of the one 
commonly used in the solution of the EW equation [13], [28]; by deflnition, it is symmetric 
under t ^ tw Hereafter we take t >t^, dropping the absolute value in the exponentials. 

One can easily obtain several limits that show the complexity of the aging behavior. 
Let flrst consider the stationary limit, which is reached when uk'^t^ ^ 1 for all n in 
(fTOl) . This condition is fulfllled when t^ ^ ti, where 

tL^T^ (11) 

is a characteristic time marking the onset of flnite size equilibration, and which will play 
an important role in the expressions below. Then, in the stationary limit, the structure 
factor becomes 

hm (5„)(t,tJ = ^(l-e-^'="^*). (12) 

If one subsequently takes the small wave- vector limit for a given At , i.e. 
I'k'^At <^ 1, the structure factor reaches the generic asymptote (Sn) = 2TAt/7. On 
the other hand, if one takes the large wave- vector limit for a given At, i.e. I'k'^At ^ 1, 
one flnds the power-law decay of the structure factor 

2T 

lim (5„)(t,t^) = — -, (13) 

uk^ Af>l 

typically found in harmonic processes. Thus, this is the typical stationary solution of the 
structure factor, displaying the At-dependent saturation in the small wave-vector limit 
and the l/k"^ behavior in the long wave- vector limit corresponding to the modes already 
equilibrated at the temperature T [13]. When At ^ t^, all the modes are equihbrated 
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at the working temperature and the structure factor shows the power-law behavior in 
f|T3|l for all kn, without the saturation regime. This behaviour indicates that the large 
wave- vectors equilibrate first at the working temperature, and then, for increasing time- 
delay, the number of equilibrated modes increase. This is obviously related to a growing 
correlation length, which will be analyzed in Sec. [3], together with its aging behaviour. 

Now, let see how the waiting-time dependence modifies the behaviour of the 
structure factor. First one notices that the small wave-vector saturation asymptote 
does not depend on the waiting time, i.e. the saturation value 

2T 

lim (S„)(t,t^) = — At, (14) 

is ti„-independent but time-delay-dependent. On the other hand, the waiting-time is 
most prominent in the asymptotic power-law regime, uk'^At ^ 1, where one finds 

lim {Sn){t,t^) = ^e~^'"'"'- + ^(2-e-^'"'"'A. (15) 

From this expression one observes that the asymptotic power-law regime in the large 
wave-vector limit changes from the (Tq + T)/ {'juk'^) asymptote at small waiting-time, 
to the 2T/ (72^^^) asymptote in the long waiting-time limit. For a fixed time-delay value 
such that At > t^, and for modes such that i/fc^At ^ 1, the power-law behaviour 
changes from (To + T)/(7z/fc^) at lower values of k^ to 2T/(7z/fc^) at higher values, 
indicating that the large wave-vectors vk'^t^, ^ 1 are equilibrated at the working 
temperature, while the other modes still reflect the initial condition. This also implies 
the existence of a waiting-time-dependent mode, k^ ~ (z/t^)~^/^, which separates these 
two regimes, allowing us to write the structure factor, in the asymptotic time-delay 
limit, in a scaling form: 

lim {Sn){t,t^) = k-'' S{k/k^). (16) 

In Figs. [2111] we display the analytic results explained above. In Fig. [2] we show 
the T = 1 evolution of the structure factor {Sn) after heating a perfectly fiat initial 
condition, i.e. from Tq = 0. The first panel illustrates the time-delay saturation value 
and the power-law decay at large where the change between the two asymptotes, 
which are proportional to Tq + T and 2T, is clear. Figure [2] (b) shows the crossover 
at ku, ~ (z/tu;)~^/^, between the two asymptotes using different waiting-times and the 
large time-delay hmit At ^ ti. In this case, i.e. T > Tq, the equilibrated asymptote, 
proportional to 2T, is the upper dashed line. In Fig. [3] (a) we show the dynamics after a 
quench from Tq > T; the approach to the saturation value at small kn is non-monotonic. 
In the asymptotic time-delay limit shown in panel (b), (Sn) crosses over from a higher 
value asymptote at small wave- vector to a lower value at large wave- vector, proportional 
to 2T and corresponding to equilibration, at the waiting-time dependent value kw in 
equation (flGll . In Fig. |l]we study the dependence on the initial condition, Tq, finding 
that the height and width of the bump increase with Tq. 
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Figure 2. T = 1 evolution of the structure factor from a flat initial condition 
(To = < T). (a) tu, = 10 and the curves correspond to the time-delay values 
At = 1, 10, 10^, 10^, 10'', and 10^ from bottom to top. The upper and lower 
dashed lines are 2T/{'yvk^) and T/(7z/fc^), respectively. The crossover between these 
asymptotes at ^ t^, (kn ^ k^) is clear, (b) Asymptotic time-delay limit, At ^ II, 
for tw = 1, 10, 10^, 10'^, 10^, and 10^, from bottom to top. Dashed lines as in (a). 




Figure 3. T = 1 evolution of the structure factor from an equilibrium initial 
condition at Tq = 5 > T. (a) tw = 10 and different time-delay values At = 
1, 10, 10^, 103, 10^, and 10^, from bottom to top. The dashed lines are 2T/{juk'^) 
and (To + T)/ {jiyk'^). (b) Asymptotic time-delay limit. At ^ t^, for finite waiting- 
times tu, = 1, 10, 10^, 10^, 10"', and 10^, from right to left. Dashed lines as in (a). 




lo' 10'^ 10"' 10° 

k 



Figure 4. T = 1 evolution of the structure factor for different initial conditions, 
To = 0, 0.2, 0.5, 1. 2. 5, 10, from bottom to top. t,,, = 10 and At = lO''. The dashed 
lines are T/{'yi^k^) (middle) and {T -\- To)/{'ji^k^) with To = (bottom) and To = 10 
(top). 
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10° lo' 10" lo' 10* lo' lo' lo" lo' 10" 10"' 10* 10^ lo' 

At At 

Figure 5. Aging of tlie two-times correlation lengtli l{At,tw) for L = 1000. (a) 
To = 1 < T = 5 and (b) Tq = 5 > T = 1. The continuous curves are the full solution 
with tuj = 1, 10, 10^, 10'^ from left to right in (a) and from bottom to top in (b). The 
thick curves correspond to limt^^tL ' ^-^d limt^<gf^ /; the former is located below for 
T > Tq and above for T < Tq. The dashed straight lines correspond to the infinite size 
limit i — > oo taken at the outset and describe the At-^^^ growth law. 



3. Two-times correlation length 

A practical definition of a growing two-times correlation length, l{At,tyj), is given by 
the fact that it marks the transition between the saturation regime at small kn and the 
behavior at large k^- 

{Sn=l){At,t^) = lim {Sn=L/l)i^t,t^). (17) 

This is similar to the numerical study in [13] for disordered elastic lines that, however, 
focuses only on a fiat initial condition and does not take into account the waiting-time 
dependence. In that case, for t^, = and Tq = 0, a crossover from a power law growth 
to a logarithmic growth was found, which is related to the disordered energy landscape. 
In our case, without disorder, we found a waiting time crossover between two regimes 
marking the separation between non- equilibrated and equilibrated modes, as detailed in 
the following subsections. 

3.1. Stationary limit. 

In the stationary limit t^ ^ ti one finds 

lim l{At,t^) = L\/\ - e-^*/2*^. (18) 

This expression varies from limAt<t^ limt^^j^ /(At, t^„) = V^vrVAt in the growing 
regime to limAt>t^ lim^^^i^ Z(At, t^^) = L in the saturation regime. These limits 
correspond to the increases of the saturation value of the structure factor and the limit 
in which all the modes are equilibrated at the working temperature, respectively. The 
saturation time in the stationary regime corresponds to the time-delay when these limits 
are equal, At^^^j^ = 2tL, and it is independent of T and Tq. 
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3.2. Decay of the initial condition. 

Setting tw = (or, more generally, t^ ^tLj one finds 

/To (1 - e-^*/2tL)2 + T (1 - e-^*/*i) , , 

lim /(At,U =lW^A ^ ^ (19) 

which varies between limAt<t^ lim^^^^i^ /(At, t^,) = ^JSn'^v/S.tT / [Tq + T) in the growing 
regime and limAi>i^ limj^^j^ /(At, t^^,) = L in the saturation regime. Equating these 
limits one finds the saturation time associated to the growing correlation length 
At*^^Q = (Tq/T + 1) ti, which depends on the initial condition. 



3.3. Aging scaling. 

In between the initial and the late stationary growths the two-times correlation length 
ages; i.e. it depends also on t^. Equation ( |T71) can be recast in a way that makes the 
scaling solution apparent 



ai^.^'-^]- (20) 



Indeed, this equation has a unique solution for /^/t^ for each set of parameters in the 
right-hand-side. In the growing regime. At ^ ti, for Tq > T the two-times length / 
moves from the upper asymptote, that corresponds to the t^ = form in equation (IT^ . 
to the lower one, that corresponds to the t^ ^ t^ form in equation f|T8|) . at a At that 
increases with tw For Tq < T the trend reverses: / moves from the lower asymptote to 
the upper one, representing the tw <t^tL and t^,^ tL limits, respectively. The crossover 
between the two asymptotes occurs at a tit,-dependent At. More precisely, in the growing 
regime one has 

/^(At, t., T, To) = I '^^^^1^^^ ^ ^ \l (21) 

with Coo = 47r^z/. It is clear that the relative value of the prefactors depends on Tq > T 
or To < T. The prefactor cq = 2coo/ (Tq/T + 1) vanishes at T/Td <^ 1, equals Coo at 
Tq = T and approaches 2coo at T/Tq ^ 1. The behaviour of the two-times dependent 
correlation length is summarized in Fig. [5] which displays the waiting-time dependent 
growth and subsequent saturation regime for two initial conditions (a) Tq < T and (b) 
To>T. 

One should notice that / increases monotonically with At and t^j before reaching 
saturation for both Tq > T and Tq < T. This increase displays a multiplicative aging 
scahng behaviour between two asymptotes with opposite trend depending on To > T or 
To < T. These kind of results where obtained, with additive aging scaling, in the out- 
of-equilibrium relaxation of mixtures of soft spheres and Lennard- Jones particles [18] 
and the 3d EA spin-glass [21], [19] after a quench from a high temperature (although the 
finite size saturation was not reached in these numerical studies). The heating case was 
not considered in these models. 
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4. Two-times roughness 

The statics and equilibrium dynamics of elastic manifolds is usually understood in terms 
of the time, temperature and system size dependence of their averaged roughness or 
width [l]. In most studies of interface dynamics one compares the time- dependent 
configuration to the initial state, typically taken to be perfectly flat (equilibrium at zero 
temperature). The thermal averaged (one-time dependent) roughness is then 

{w^) (t) = L-'' J d'^z { [x{z, t) - x{z, 0)]') , (22) 

where x is the height of the surface and z is the position in the d- dimensional substrate 
typically with cubic geometry and linear size L. 

The initial grow of the roughness is linear with time, 

hm (w^)(t) = 2Tt, (23) 

which corresponds to a normal diffusion regime in which the beads on the line are still 
uncorrelated. Note that there is no ballistic regime since the EW equation represents 
overdamped Langevin dynamics. However, these inertial effects could be relevant in 
polymer - molecular - dynamics studies. The normal diffusion regime can be considered 
as a transient before the correlated dynamics of the interface is reached. In the single 
particle regime the roughness does not age. In the following we concentrate on the 
correlated aging dynamics of the line. 

The thermal (and disorder averaged if random interactions are present) roughness 
follows the Family- Vicsek scaling [29], which means that it crosses over from growth to 
saturation at t^ ~ [Ij: 

{w'){t)r^L^f{t/t^), (24) 

where the scaling function obeys f{y) ~ for ?/ <C 1 and /(y) ~ 1 for y 1, with 
C, P and z = (/P the roughness, growth and dynamic exponents, respectively. For 
the Edwards- Wilkinson (EW) line in 1 + c? dimensions ( = 2 — d, (3 = 1 — d/2 and 
z = 2. In the presence of disorder, ( is expected to take a 'thermal' value, (th for 
L < Lc{T) and a larger 'disorder' dominated value, (dis for L > Lc(T), both exponents 
being T-independent [1]. The other exponents, (3 and z may depend on T [H] or even 
logarithmic time-dependencies may exist p!3]. 

We wish to take into account the waiting-time t^ and consider also more general 
initial conditions. To this end we generalize the deflnition in equation (!22|) to 

I rL 2 °° 

{w^)(t,t^) = - dz {[6x{z,t)-6x{z,t^)f) = -J2{Sn){t,t^) (25) 

n=l 

with Sx{z, t) deflned in equation ([3]) and especialized to d = 1. Note that the zero mode 
is not included in the sum. 

In general, one expects the dynamics to become stationary after an equilibration 
time ti, the generalized thermal averaged roughness should then scale as in equation 
with t replaced by At. For not too short L, ti may become very long and the dynamics 
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might remain non-stationary with {w"^) depending on t^ exphcitly for t^ < t^. In [O] 
we conjectured that the scahng of the roughness in the non- equilibrium relaxation of 
infinitely long elastic lines with or without quenched disordered potentials follows the 
law 

m 



{lV^){At,t^) f{t^) T 



(26) 



with lit) a growing length (dimensions are restored by prefactors that we omit) and 
T a scaling function. For each waiting-time this form approaches a stationary growth 
regime {u?) ~ £''(At) when t^; ^ At ^ if T{y) ~ y'' for ?/ ^ 1. It is also reasonable 
to assume T\l{t) j l{tw)\ ~ i"'(At) for At ^ t^, that leads to a stationary growth of the 
averaged roughness at very short time-delays. This result is explicitly realized in the 
power law case. One may extend this conjecture to apply to manifolds with internal 
dimension D in a transverse space with dimension d. The functional form of the growing 
length, the scaling function, JF, and the values of the exponent are expected to vary 
from case to case. 

In a series of numerical studies one established that, in the numerically accessible 
times, the quenched dynamics of a disordered 1 + 1 lattice model [U [9] satisfies the 
scaling in equation (1261) with t{t) ~ t"/^ and a rather sma// exponent. A crossover 
to a logarithmic time-dependence [13] is not excluded although it was not seen in the 
out-of-equilibrium relaxation. It was shown that the roughness ages, by crossing over 
between two asymptotes, for At ^ t^ and At <^t^, thus having 

(w^) ~ ci,2(T) At"(^) , with a{T) < (3ew = 0.5 , (27) 

and different proportionality constants. The waiting-time dependence appears in the 
way these asymptotes are approached. aiT) is a generalization of the growth exponent, 
P, in surface growth literature, and a(T) < Pew, with Pew the Edwards- Wilkinson 
value, reflects that quenched disorder slows down the dynamics. 

In order to better understand this aging behaviour, trying to separate the effects due 
to disorder from those related to the intrinsic elastic character of the line, we present 
here results for the EW case in 1 + 1 dimensions without disorder. In this case the 
two-times averaged roughness is given by 

{w^){At,tJ = ^ Yhn{At,t^) + ^ Va„(At,t^), (28) 

n=l n=l 

n^a„(At,tJ = 2 (l - e'"^^) (l - e^^t) + (l - e""^^) e^^t, (29) 
n%^{At,t^) = 1-e " e " (30) 



where we used the definitions 

w;2 = ToL/(127z/), = T L j {\2^v) . (31) 

The averaged two-times roughness can be written in terms of the scaled times 



2\/A ^ _ / 2\ ( tv\ _ 2\ f 



{w'){At,t^) = {w') (^^,^J = K) [j^^-^) (32) 
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while temperatures appear separately, through and w^. We first focus on the time- 
dependence of the roughness of infinite lines, L —>■ oo, before testing the dynamics at 
different asymptotic limits, and we later reverse the order of limits by considering finite 
lines, L < oo. We display the analytic results in Figs. [MHlin which we approximate the 
series in the analytic expressions by using finite sums, Yl'^=o ~^ J2n=o^ with M = 1000 
in all cases, except for lines with L = 3000 for which we used M = 500. 



4.1. Infinite system size. 

By taking L oo at the outset ti diverges and 



lim {w'){At,t^) = J^ 



(33) 



which can be written in the scaling form 

lim {w') (At, t^)=tl/'w'(^), (34) 

and admits the scaling form in equation fl2^ with i(t) ~ t^/^. Comparing now At to 
a long waiting-time, that is to say taking At ^ t^, one finds that the waiting-time 
dependence determines the crossover between two square-root dependencies in At with 
different, temperature-dependent, prefactors: 

{w )(At, t J - j ^^.^^^ (35) 
with cf = [T + To{V2 - l)]y^2/{^n^ and = 2T/^/^^. Note that cf < d£ for 

2 2 

T > To, the two constants are identical at T = Tq, and Cq > for T <Tq. The two 
trends are shown in Fig. [7] (b). These results are of the generic form proposed in [9], 
see equation with a = (3ew = 1/2 independently of temperature, i{t) ~ t^/^, and 
temperature dependencies of the constants made explicit. 



4-2. Finite lines 

For finite L, by considering the long waiting-time limit t^ ^ t^ the stationary Family- 
Vicsek scaling form (l24l) is recovered. Indeed, one has that 



,Um(..^)(AM.,) = .iEAj(l--"'^''='')- (36) 

n=l 

This form for the stationary roughness contains the growing regime where 
limAf<f^ limf^^^t^ (w^) = cJ^^(T)At-'^/^ as in eq. ([35]), and later crosses over to the 
saturation value 

TL 

lim lim (^i;2)(At,t^) = 2wl = (37) 

On the other hand, the averaged roughness at tt„ = for finite L, or more precisely 
for t^ -C ti,, also follows the familiar Family- Vicsek scaling form fl21|) . In this case, the 




10' 



10 

At 



10 

At 



Figure 6. (a) Stationary roughness (t^ S> t^) as a function of At. Thin lines 
correspond to different system sizes L = 100, 300, 1000, 3000, from bottom to top all 
at T = 1. Thick lines correspond to L = 1000 and T = 2 (bottom) and T = 5 (top). 
Note that the saturation time does not depend on T. The dashed line is the limit 
L — + CX3 for T = 1. (b) Roughness with = 0, T ~ 1, and different initial conditions 
corresponding, from bottom to top, to Tg = 0, 1, and 5. In the t^, ^ limit, the 
roughness goes to its stationary solution, also given by the curve in the middle. Dashed 
lines correspond to limi^oo(w^) given by equation (|35p. 



roughness crosses over from a growing regime limAt<t^ lim^^^^^ = Cq ^(T, Tq) At^/^ to a 
saturation value which depends on the initial condition, 

hm hm {w'){At,t^)=wl + wl = + (38) 

Note that the saturation value reached in the stationary regime, after having taken 
tw ^ thi is twice the saturation value obtained from the flat initial condition Tq = 0. 
This result demonstrates how important it is to be careful with the choice of to ensure 
that one has reached the stationary regime when using numerical simulations. 

When considering all the waiting-time dependence the roughness interpolates 
between the different limiting values given above. For sufficiently large but finite systems 
it is possible to find a well defined aging of the growing regime when both At and 
are smaller than the saturation time. In this case one recovers the scaling behaviour 
found in the infinite size system, equations (133!) and (l35l) . For finite L though this 
scaling form terminates at a characteristic time-delay, tx{tw) oc ti - see below - and 
{w"^) later saturates at a waiting-time-dependent value, which interpolates between the 
two limiting saturation values, equations ( l37l) and ( l38l) . These results are the same 
as the ones obtained using the reverse order of time limits, see flTTI) and fl55]l . The 
t^-dependence in the saturation regime is shown in Fig. [HI 

Although in the previous analysis we used tL as a reference time to determine 
the long time-delay limit. One can be more precise and define a waiting-time 
dependent saturation time, txitu,)- Equating the growing and saturation behaviour 
of the roughness, one can extract the saturation times at two extreme values of t^,'- 

' L\ (39) 



T + Tq 
T + To(v^ 
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Figure 7. Aging of the two-times rougliness, = 0, 10°, 10\ 10^ 10^, 10'', 10^. (a) 
L — 1000, To = 5 and T — \, increases from top to bottom. The upper and lower 
dashed lines correspond to L ^ cx3 with tw = Q and tw ^ thi respectively, (b) L ^ oo; 

increases from top to bottom in the lower set of curves (Tq ~ 5 and T — 1) and 
from bottom to top in the upper set of curves (To = 1 and T = 5). 

Io'l — — — — !"i — a 




10° lo' 10" 10"' 10^ lo' lo' 

Figure 8. Aging of the roughness at saturation, limAt>ti w^(At, tj„). T — 1 and 
To — 0, 0.2, 1, 2, 5 from bottom to top. Dashed lines correspond to {T + To)L/ {12 jiy). 

t.{t^:^tL) = ^L\ (40) 

Both results are proportional to t^, and txitw) interpolates between these two values. 
This justifies the use of as the reference time-delay before saturation; the prefactors 
are finite for all T and Tq, the former is larger (smaller) for T < Tq [T > Tq) and become 
identical at T = Tq. The saturation time and its dependence on L are visible in Figs. [6] 
and [7| (a). 

4-3. Scaling relations 

Here we present the main results concerning the aging properties of the roughness. We 
have shown for the aging regime that 

tx ~ a{T) L^, z = 2, 

Mak<*. ~ c(T,To,AVUt^, /3 = l/2, 
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see equations (15^ and (HUl) , equation (pTj) and equation , respectively. The prefactor 
c{T,To,At/tJ approaches c^{T){At/t^Y/^ and c^'(T, To)(At/t^)i/^ for At < and 
At tw, respectively. These relations imply 

- aiT)iwly/< ^ [aiTfWJ'^^ , (41) 

where we have used that {w'^)At^tx ~ w"^ for simplicity. From matching the end of the 
growth regime with saturation at t^ <^ At = t^ one has 

c(T,To,At/tJtf (T,To)tf and P = C/z, (42) 

(the latter condition is satisfied by the values of the exponents found). Thus one finally 
has that 

a(T)~f/^ = cf (T,ro). (43) 

Note that the exponents take simple values in the EW case but, in general, they can be 
T-dependent (not () [TTj . 



5. Mean-squared and center-of-mass displacements 



The two-times roughness, equation ( l25l) . may also be written as 

{w^){t,t^) = {B){t,t^) - {D){t,t^) 

= (^[xiz,t) - x{z,t^)f'j - - xit^)]^^Y (44) 

where (B) and (D) represent the averaged mean-squared-displacement of the differential 
line segments and the center of mass of the line, respectively, see Fig. [H This relation 
simply states that the roughness is a measure of the fluctuations around the center of 
mass of the line. 

It is simple to show that the center of mass diffuses normally. Indeed, integrating 
equation ([1]) over the line length one has dtx{t) = ^'{t), with ^'{t) = dz ^{z,t), 
(m) = 0, {em'{t')) = '^6it-t%and 

{D)it,t^) = {D){At) = '^At. (45) 

The diffusion constant is an inverse function of the line length. 

In previous studies of the elastic line out-of-equilibrium dynamics [HI [16] one focused 
on the mean-squared-displacement that, in the case of the EW line, is just given by 
{B){t,tw) = {w'^)(t,tw) + 2T/(7L) At. At short time-delay and B are practically 
identical while in the saturation regime the displacement is just given by the normal 
diffusion law. 
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Figure 9. The wave-vector dependent correlation function defined in equation (|47|) 
and given in terms of {w^) by equation (l48l) for a Gaussian process. {Cq) at Tq = 5, 
T = 1, g = 0.1 and tyj = 0, 1, 10, 10^, 10^, 10^ and 10^. (a) L = 1000. The thick 
dashed lines correspond to the limits L oo with — (lower curve) and 3> 
(upper curve), (b) L — > oo limit; the rest of the parameters are as in panel (a). The 
inset shows the same curves in a different scale, i.e. In(Cg) vs. Ai^/^. 



6. The incoherent scattering function 

The dynamics of glassy systems is usually analyzed in terms of the wave-vector 
dependent incoherent scattering function: 

N 

{Cg){t, t^) = N-^J2( e'^>^>W-"^>(*-)] ) (46) 

i=l 

with N the total number of particles, ri(t) the time-dependent position of particle i 
and q a wave-vector. Cq is measured numerically and experimentally. In the context of 
elastic lines, one defines 

(C,)(At,t^) = L-^ j dz ( e^9[^^("'*)-'5^(^'*-)] ). (47) 

In the EW case the displacement in the exponential is a Gaussian random variable and 

(Cg)(At,t^) = e-^<-')(^*'*-). (48) 

The incoherent scattering function (Cq) is simply related to the roughness, (ti'^), 
analyzed in Sect. HI Figure [9] displays the time-delay decay of {Cq) at g = 0.1 using 
several waiting-times. For finite lines the saturation in {uP') is attained at sufficiently 
long At and the correlation reaches a waiting-time dependent plateau with its height 
increasing with t^, Figj9] (a). For sufficiently long lines the saturation regime can be 
pushed beyond the observed time-delay window as exemplified in the limit L —>■ oo, 
Fig. [9] (b). In Fig. [10] we display the dependence of {Cq) on At for fixed and different 
values of q. One notices that the small q correlations saturate in the At window while 
large q correlations relax to zero. This is to be expected since g^w^ scales as q^L and 
the effect of decreasing q is like decreasing L. Note the similarity between these plots 
and light-scattering measurements in clay suspensions (laponite) [17]. 
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Figure 10. {Cq) for L = 1000, = 10^ and different wave-vectors q as indicated. 
To = 5 and T = 1. 



7. Response functions and FDT 

We now compute the linear response function of several observables related to the two- 
times correlations studied above. 



7.1. Center of mass response 

In order to evaluate the linear response we have to switch on a perturbing field coupled 
to the observable of interest. Let us start with the linear response function associated 
with the center of mass diffusion, (x^) . The effect of a perturbing field, h, apphed on 
the center of mass after time ty^ is described by the term 

n'^ = -h Lx{t) 9{t - t^) = -h I dz x{z, t) 9{At) (49) 

Jo 

that is added to the energy. Calling x^{t) and x{t) the center of mass position with and 
without field, respectively, the linear response function is 

{X''){t,tj = ^{x\t)-x{t)) (50) 
which depends only on the time difference and satisfies the FDT for any t and ty,, 

(L>)(At) = 2r(x^)(At). (51) 

7.2. Roughness response 

The energy contribution of a field conjugated to the roughness is 

n""'^-h / dz [x{z,t) -x{t)]s{z)e{At). (52) 
Jo 

As usual, s{z) are i.i.d. quenched random variables taking values s{z) — ±1 with equal 
probability: {s{z)) — and {s{z)s{z')) — S{z — z'). The associated response function is 

{X'"''){t,t^) = ^ Q^'^dz [dx'^{z,t)-dx{z,t)] siz) 
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w^(At,t ) w^(At,t ) 



Figure 11. Violation of the FDT in the aging regime of the EW equation. The 
parametric plots = l/(27e//) are constructed from the scaled variables defined 
by w"^ = tw^^^{w'^) and x™^ = tw^^'^ {x^') ■ (a) The parametric plot in linear scale 
showing the departure from the FDT (dotted line). Upper and lower sets of curves 
correspond to Tg = 1, T = 5 and Tg = 5, T = 1, respectively. Thick dashed lines 
represent the L oo limit. Thin lines are for L = 1000 and different waiting- 
times, and show the finite size signature in the FDT parametric plot. The inset 
shows the large scale violation of FDT with two straight lines, which corresponds to 
effective temperatures larger (for Tq > T) and smaller (for Tq < T) than the working 
temperature T. The initial FDT regime is not clearly observed in this scale, (b) The 
parametric plot in log-log scale showing that the violation of the FDT is given by 
equation ([63|l at At:$>tyj. 



2 °° 



= lY.i[^n{t)-Cn{t)]Sn), (53) 
n=l 

with Sn defined through 6s{z) = s{z) —s = Yl^=-oo^n^^^"^ ■ Here (. . .) also indicates 
the average over the s„ distribution. One finds 

(x"')(At) = fl -e-'^'^'^*) (54) 



where we used (s^) = 1/L. Interestingly enough, the linear response is stationary for 
all tw while the roughness is not. Therefore, the FDT is not respected for waiting-times 
tw ^ ^L, and its modification is discussed in Sect. 17.41 and Fig. [TTl In the stationary 
limit ty^ ^ ti the roughness becomes stationary and the FDT holds, i.e. 

hm {w'){At,t^) = 2T(x"')(At). (55) 

This statement also implies that the FDT does not hold for t^ = and T = 0, pointing 
again that one should be careful with the choice of t^ and the stationary limit. 

7.3. Mean- squared- displacement response 

The effect of a perturbing field conjugated to the mean-squared-displacement is 
represented by 

n'^ = -h [ dz x{z,t)s{z) 9{At), (56) 
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Figure 12. The effective temperature Te// clraracterizing the violation of the FDT in 
the EW equation, (a) Evolution with the rescaled time At/t^ showing that FDT holds 
when At <C and it does not hold when At 3> . Upper and lower curves correspond 
to To = 5, r = 1 and Tq = 1, T = 5, respectively. The dashed lines represent the limit 
At > tni in equation ([63]). (b) Teff{At > t^) vs. T for To = 0.3 (continuous 
line) compared to the data in Fig. 25 in [1^ describing the effective temperature of 
independent elastic lines (with the Josephson correction to the elasticity) moving in a 
quenched random environment. The inset shows j/At>i„ = liiiiAt>t,„ T'/Te//(T) as a 
function of temperature for Tq = 1, T = 5 (lower curve) and Tq = 5, T = 1 (upper 
curve). The dashed line is the limit T ^ Tq. 



with the random s{z) distributed as above. The hnear response function is defined as 

= ^(|^ dz [x''iz,t)-x{z,t)]s{z)'j 

= {x''')it,t^)+^-^[x\t)-xit)]. (57) 

The last term of this expression represents the center-of-mass response to a field of 
intensity h' = h/(s). Thus, in the long waiting-time limit, {x^) is also stationary and 
simply related to the roughness and center-of-mass responses, 

(x^)(At) = (x-')(At) + (x^)(At). (58) 

In this case, the FDT is not satisfied in general but in the stationary regime it is: 

lim (5)(t,t^) = (fi)(At) = 2T(x^)(At). (59) 

7.4- FDT and effective temperature 

Since FDT holds for t^ ^ ti, it is interesting to study the violation of the FDT at finite 
t^. To this end we use the L ^ cxd limit, where the roughness takes the scaling form in 
equation (l33!l . and 



lim {x-'){At) = t]l-r with = (60) 

L^oo \tu, \tw J V vr7V t^ 
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Figure 13. The effective temperature Tg// as a function of time-delay for tw = 
1, 10, 10^, 10^, 10"* and 10^ from top to bottom; L = 1000, Tq = 5 and T = 1. The 
dashed Hne is linit^^At lin:iL-.oo 7e/y [the same dashed curve as in Fig. [T^ (a)] and 
justifies the shoulder for small t^- The inset shows Teff{At ~ 10^ 3> t^jtw) as a 
function of t^', it is clear that at sufficiently long one recovers T^ff = T. 



Once we eliminated the ti/^ diffusive factor one can associate the ratio between w"^ and 
X'^' with an effective temperature [8l [TO]. 

= T^w'. (61) 



which, in this case, depends on T, Tq 

T-To 



At and t^- From equation fl33l) we obtain 



e// 



1 + 




2At ^ 



1 + 



At 



1 + 



2t 
At 



(62) 



In the case T = Tq one recovers Tejj = T as expected. One can check that Tgjj > T or 
Te// < T whenever Tq > T or Tq < T. In the two extreme cases At <^ tw and At ^ t^, 
one finds two waiting-time independent values of Tf.ff: 

T At < t^. 



T+(To-T) 1 



At > t. 



(63) 



see Figs. [TT] and [12], indicating that fast modes are equilibrated while the slow ones 
are not. These results are very similar to what has been found numerically for the 
dynamics of elastic lines in a quenched random potential [8] and in models of interacting 
elastic lines in quenched random environments that describe the vortex glass in high Tc 
superconductors [10], see Fig. [121 Note that the waiting-time dependence in T^jf only 
marks the crossover between the two asymptotic regimes. At still longer At such that 
At 3> ta; ~ ti, the FDT result, T^ff = T, is recovered. This is shown in Fig. [131 a result 
that, once again, resembles what was found in laponite [30] . 



8. Fluctuations 



Upto this point we studied the scaling properties of several two-times quantities averaged 
over the length of the line and thermal noise. A more refined investigation of interface 
dynamics, currently done theoretically and experimentally, deals with the fluctuating 
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quantities such as the distribution functions of the width of heights [20], the density of 
local maxima or minima of heights [311 132] . the statistics of first passage times, etc. It 
is also clear by now that to gain a complete understanding of glassiness one should also 
understand the dynamic fluctuations [211 122] . In the glassy context one expects that 
the competition between different time or length scales in the system reflects in the way 
the probability distribution functions (pdfs) behave. 

Racz proposed that, for elastic systems, the scaling form of the distribution function 
characterizing the roughness fluctuations may serve to classify the systems into different 
universality classes [20] . For a given system size the pdf of the roughness, Pl(w^), scales 
as [33| 

wIPl{w') = $ (^) , (64) 

in the saturation regime, and the form of the scaling function was shown to be useful 
to differentiate between the EW and KPZ universality classes [20]. The full stationary 
pdf of the EW roughness was also computed, showing the same scaling behavior as in 
equation ( IMIl . These works considered only the t^ = case with a flat initial condition 
(T = 0). Recently, a simulation study of a disordered elastic-line model defined on the 
lattice demonstrated that the scaling form of the distribution function must be modified 
to include the aging effects [5] • 

In the body of this Section we analyze the thermal noise-induced dynamic 
fluctuations of the two-times quantities defined previously during the aging relaxation. 

8.1. Roughness distribution 

Rewriting the roughness in terms of u{z,t,tw) = x{z,t) — x{z,tw), the pdf of the 
two-times roughness, Pl^w"^), is given by 

Pl{w^) = J V[x]V[x'] 6 _ (u^^t, t^) - u{t, t^)')] p (x, t; x' , t^) , 

where V[x] is the measure over x configurations; x and x' represent the configurations 
at time t and t^, respectively; and p {x,t; x' ,tui) is their joint probability density. The 
Laplace transform Gl{X) = da Pl{o:) e~^°, can be written as the path integral 

Gl{X) = [ V[x]V[x'] p {x, t- x', t^) e-^("'(*'*»)-^^(*^') . (65) 



Using the independent Fourier modes defined in equation ([3]) 

/oo 
\[dcndcldd^ddlp' [c„(t),c„(t^)|c„(0)] e-2^l^"W~^"(*-)l', (66) 
n=l 

where A/" is a normalization factor ensuring G'l(O) = 1 at all times. The quantity 
p [c„(t), c„(tiu)|c„(0)] is the joint probability density of having c„ at time t and = 
Cnitw) at time t^,, given the initial condition c„(0). This quantity can be expressed 
as p[c„(t),c„(t^)|c„(0)] = p[cn{t)\cn{ty,)]p[cn{tyj)\cn{Q)], where p[cn{t)\cn{t')] is the 
conditional probability of evolving from c„(t') to c„(t) in the time interval t — t'. This 
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states that the joint probabihty is simply the product of the probabihties of evolving 
from the initial condition to the intermediate configuration c„(t^), and from there to 
the configuration c„(t), and satisfies 

P[c„(t)|c„(0)] = j dc'„dc'lp[Cn{t),cM\Cnm- (67) 

The conditional probability p [c„(t)|c„(t')] is given by the complex Gaussian 



p[Cn{t)\Cn{t')] = ^^^2lt_t,^ ^ ' (6^ 

where 



|c„(t)-c„(t')e-''4('-'')|^ 



al{t-t') = -^\l-e-'^'^'^('-''^]. (69) 
Therefore, the joint probability function for the Fourier modes at times t and t^ becomes 

I i2At|2 I ,2. |2 

\c„(t)^cn(tw)e-^''n^-t\ c„(tu,)-c„(0)e-''''>i*™ 



e 2^J^At) g 2<T|l(t„) 

P [C„(t), C„(t^)|c„(0)] = /9^N2^2rAfW2rf ^ • C^O) 

(2vr)^<(At)(T^(t^) 

After some algebra one finds the normalization factor Af = Y[^=i 167r^cr^(At)cr^(t^) and 



OO 



g l + Au)^a„(At,tu,) 



„=i 1 + Aw;4a„(At,t^) 

with the coefficients a„ defined in equation (129!) . This is the two-times generalization 
of the result in ^23j, including arbitrary initial conditions. The averaged two-times 
roughness follows from Gl as {w^){At,tw) = —dxGL{X,t,t^)\^^Q, which allows to 
recover the result in equation ( l28l) for the initial condition |c„(0)p = To/^Ljukn). 
Since the two-times roughness pdf is given by 

/ioo j\ 
— e^-'Gi(A,t,tJ, (72) 

we can extract its scaling properties from the ones of G'l(A) in equation ( |7T1) . Using 
y = Xw"^ and 2|c„(0)p/u7^ = 6iyo/(vr^n^w^) = 6/(7r^n^) Sg we find 

At 



«^L^L(^^) = $(f^;=^,r'^o), (73) 



!/SQbn( 



At tw ' 



f \f f \ ri'^ rill ^ p i+'"'"(#'^f) 



with the coefficients a„ and 6„ defined in equations fl29|) and fl30|) . For a fiat initial 
condition, Tq = and Sg = Wq/w"^ = Tq/T = 0. With a very similar calculation to the 
one explained in [23] for the stationary case, we rewrite the function $ as 

\ tL ti J an an -am 
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which is essentially the same result in [23] but with two-times dependent coefficients 

{At/tL,t^/tL). 

By using now a different independent variable, x' = w"^/ (w"^), one has 

{w')Pl{w') = <^'(^;^,'-^,sI], (76) 



I 2,1 ( 



At tw ' 
L 



' At 



( ^ ,A fr e ■-'-■(I. 

^ ' (1 • (1 • V i-ioo 27ri i i 1 + 



and 



/ '^n ,/ bn (7Q') 

Equations (I73|) and (17^ are the generalization of equation (15) in [23] that takes into 
account the aging regime. Equation ( 1761) is a rewriting of the latter using the more 
convenient normalized variable w'^/{w'^). The parameters are, in both cases, At/t^, 
tw/tL and Sg. In the growth and aging regime in which At/tyj is finite and the two 
parameters are very small, i.e. At/ti, tw/tL <^ 1, one formally has 

^ = e, w' = G{e\ {w^)=t]l^g{e). (79) 

We can now easily exchange At/ti and t^/ti by the more convenient set w"^ and 
{w^)/wl^. First, we exchange At/t^ and t^/tz, by At/t^ and t^/ti. Second, on the 
one hand At/t^ is an exclusive function of w"^. On the other hand, using the results in 
Sect. 14.31 one can show that 



ti V "^.^ 



oo 



(80) 



The factor w'^ enters the last equation, but we can ignore it by redefining the scaling 
function. We used the EW exponents but this relation can be easily rewritten for generic 
j3, z and C,. We thus have 

(^^)P,(^2) = { w\ T, To^ (81) 



as proposed in [9] for the generic disordered case. 

Let us list and illustrate in some plots different trends in the scaling function $ 
evaluated for the fiat initial condition. In the numerical evaluations we approximate the 
infinite sums and products, as in fl75]) . by finite sums and products with different cut-off 
values. Ml and M2 respectively. The main panel in Fig. [H] (a) shows the time-delay 
evolution of the scaling function $ for t^ = ^ and fixed system size. This corresponds 
essentially to the results obtained in [23], and shows how the pdf is broader for increasing 
At until the saturation regime is reached. The inset shows the system size dependence, 
indicating that at fixed t^ the function $, tends to a delta-function in the infinite size 
limit. In Fig. [14] (b) one observes how the pdf is spread at fixed L and At while increasing 
the waiting-time. These results correspond to a fiat initial condition Tq = 0. 
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w /w^ w /w^ 

Figure 14. Scaling function At/i^, f^/ii, 0) evaluated with Mi — 20 and 

Ma = 40. (a) L = 1000, = 0, and At = 10^, 3 10^, 10"*, 3 lO"*, and 10^ (from left to 
right). The inset shows a comparison between $ for different system sizes, L = 1000 
and L = 300 as indicated, at the same waiting-time i^, = 10^ (b) L = 1000, At = 3 lO'^ 
and tn, = 1, 10, 10^, 10^, and 10'' from left to right. 



Since it is hard to compute the scahng function $ for a Tq > initial condition we 
just present the skewness and kurtosis, 

a = %, .= 0-3, (82) 

respectively, with the centered moments defined as fi2 = (w^) — {w^)'^, /is = (w^) — 
3(w^) {w^) + 2(m;2)3, and /i4 = {w^) - 4(w;^) (w^) + 6{w^) {w^ -^wY■ The moments of 
the roughness pdf are given by ((w^)™) = (—1)™ d\mGL{\)\^^Q] then after some algebra 
one finds 

Figure [15] displays the time-delay evolution of the skewness and kurtosis for the cases 
T > To = [Fig. [15] (a)] and T < Tq [Fig. [15] (b)]. One can observe that the pdfs 
are broader and more asymmetric until saturation. Generally, both a and k age with 
a similar At and t^ dependencies as the correlation length / or the averaged roughness 
(w^). The peculiarities are that the asymptotes corresponding to the infinite size limits 
grow as At^/^ for the skewness and At^/^ for the kurtosis. The approach to saturation 
is non-monotonic, showing a bump around At trc{t^). Finally, Fig. [16] shows the 
To-dependence of the skewness for t^, = 10^ ^ ti (the kurtosis behaves in a similar 
way). 

8.2. Center- of -mass displacement distribution 

In Sect. [5] we showed that the line's center-of-mass undergoes Brownian motion with 
mean a;(0) and variance {D){At) = 2T/{'~fL)At. Thus, the center of mass position is 
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Figure 15. Two-times skewness a{At,tyj) and kurtosis K{At,tyj) for the distributfon 
function Pl{w^) (the kurtosis is rescaled by a factor 10 for clarity). The sums are 
truncated with M3 = 100. T = 1, L = 1000, and = 1, 10, 10^, 10^, 10^, and 10^ as 
indicated. Different initial conditions correspond to (a) Tq = and (b) Tq = 5. The 
thick lines correspond to the i^, ^ 00 limit. 




At 

Figure 16. Two-times skewness a{At,tw) for the distribution function Pl{w'^). The 

sums arc truncated with M3 = 100. T = 1, L = 1000, and To = 0, 1, 2, 3, 4, and 
5 from top to bottom. The waiting-time is tw = 10"^. For comparison, the tw = 10^ 
curve, which is the same for all initial temperatures is also included as a thick line. 



Gaussian distributed and, after a change of variables, the distribution of the mean- 
squared-displacement of the center-of-mass is 

-D/(2{D» I Zl / \ 

which can also be written in the scaled form 
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8.3. Mean-squared displacement distribution 

The mean-squared-displacement satisfies B{At,tiu) = w'^{At,tiu) + D{At). w"^ and D 
are independent variables (note that the roughness is independent of the zero mode), 
then the probability function for B can be formally obtained from the inverse Laplace 
transform of the product of Laplace transforms, 

—e'^KLiX,^t,t^), (86) 

-oo ^TTZ 

where Kl = GlJl is the Laplace transform of Pl{B), Gl{\, At,tuj) is given by 
equation (17T!) . and Jl(A, At) is the Laplace transform of Pl{D), 



Jl(X, At) = , ^ = . / . (87) 



Then one finds that 



oo 



KAX, At, tj = , TT — ^ . (88) 

^ ' v/2 [{B) - {w^)] X + 1 j^}^ l + Xwlar,{At,t^) ^ ' 

We now define y" = X{B), Aw^a„(At, t^;) = y"a'^{At,tw), and 

// I At tj^\ (3„ j At tj^\ 6„ (QQ\ 



{b)Pl{b) = ^[4,;^,'-^,sI]. (90) 



and we obtain 

XBY tL'tL 

The scaling function is given by 

J, (j'^ r'°° 'bLgi-^' n \ '"^"^^"^ (91) 

V 'tL'tL'V J-^ 2m J-iv'l + 2C>"l+!/'<(Ai,i.)' 
with 

The pdf of the mean-squared-displacement B can also be written in a scaling form 
similar to the one found for the roughness. 

8.4. The incoherent scattering function 

In order to obtain the pdf of the incoherent scattering function we need to use the 
moment-expansion: 

P{z)= j ^e-Xn) (93) 



with 

*)-E^^ -E^(^^)- (94) 



p\ duP 



—in] 
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This expression assumes that the series converges and the moments exist. The moments 
{Cl)=L-P J dzi ...dzp e"^ ^''■'=i<^> '^'^^^ (95) 



of Cq are 



where the function 

2 °° 

{H)iz;t,t^) = - Y,{Sn)iAt,t^)e'''" 

n=l 



i-L 

-1 



L / dz' {[6x{z' , t) — 6x{z' , tw)] [6x{z' — z,t) — 5x{z — tw)])(96) 



is the two-times generahzation of the height-height correlation function 

Using the fact that {H){0;t,tw) = (w^)(At,t^) and equation (jUj), one can show 
that (C^) = {Cg^Ip, where the g-dependent function Ipit.t^) is given by 

Ip{At,t^,q^L) = L~P j dzi ...dzp e'^ ^L'=i;.^.'<^>(^'-v;^*.*-). (97) 

From the moment-expansion, using u' = u{Cg), the fact that Tq^L is a function of (C^) 
and calhng T^q^L = (Cg), one can write the pdf of the incoherent scattering function 
in the scaled form 

a, At t, 



{Cg)P (Cg) = ( ^, ) , (9^ 



with 

00 



Although in the last expression the functional form of B(x) is not evident, the scaling 
properties are clear. 



8. 5. Fluctuations of the response functions 

In quadratic models as the EW elastic line the response functions do not fluctuate. 
This can be easily shown as follows. Take the roughness integrated linear response (l53l) 
without the thermal average. Replacing c^(t) and Cnif) by their functional form, and 
using the fact that the two copies evolve with the same thermal noise, one has 

00 

X"'(t, t^) = 2L ^ 4 e-'^'^" . (100) 

n=l 

This result depends on the random fields s„ but it is independent of the thermal 
noise. For fixed random fields this quantity does not fluctuate and the pdf of 
is a delta function. Similarly, one can prove that the displacement and center-of-mass 
responses are delta-distributed. The same 'trivial' result was found for the ferromagnetic 
coarsening in the O(A^) model with N ^ 00 |27j . 
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9. Summary and conclusions 

We studied the averaged and fluctuating dynamics of the Edwards- Wilkinson elastic 
line in one transverse dimension. 

Firstly, we analyzed the evolution of correlation functions in terms of the different 
time scales involved in the problem: the total time, t, the waiting time, t^, and the 
saturation time t^. In particular, for the two-times roughness we found {w"^) ~ 
J^w^ [At/tw,tuj/tL,TL,TQL] ~ J^ii,2 [t/tuj,tu,/tL,TL,ToL] (the scaling function in the 
third member is not identical to the one in the second member but we use the same 
name to simplify the notation). As mentioned earlier, the problem can also be analyzed 
in terms of associated length scales obtained from the growing correlation length 
i{t) = 47r^i/t^/^. Then, the relevant length scales are i{t), i{t^), and ^(ti) ~ L. This 
is reflected for instance in the scaling form of the structure factor in the asymptotic 
time-delay limit (flGj) . using ~ iit^)^^. The scaling form of the two-times averaged 
roughness that describes the aging, saturation and equilibrium regimes can be written 
as 

m i{t. 



{w ){At,t^) ~ J^^2 



TL,nL 



with e{t) = in^ut^/^ and ti = L'^/{47i^ u) in the 1 + 1 EW case. This form extends the 
proposal in fl26l) to include another scaling variable and thus describe the dynamics of 
flnite lines. It then generalizes the celebrated Family- Vicsek scaling [29] to include the 
preasymptotic non-equilibrium regime. The aging regime corresponds to ^(t^) <C ^(^l) 
and the function J^^2[i{t) /i{t^),0,ToL,TL] ~ fj^[i{t)/i{t^)], leading to ([26]) with all 
the temperature dependent asymptotic properties already detailed in the central part 
of the manuscript. The saturation regime is reached by taking t ^ t^ at flxed t^', this 
means (-{tw) / l{t) <^ 1 and i{t)/i(tL) ^ 1- Finally, the usual stationary equilibrium 
regime corresponds to ^ ^l, and for a power-law growth one recovers the Family- 
Vicsek scaling {w^) ~ L^T^2[i{At)/L]. 

We also showed that the two-times dependent correlation length deflned through 
the dynamics of the structure factor, satisfles a similar scaling law 

m i{t^) T' 



(102) 



.^(t«,)'^(tL)'To_ ■ 

Note that although this expression gives the full aging behaviour of l{t,tw), it can 
be completely rationalized using the simple length scale £{t). For instance, in the 
aging regime, one flnds that regions with i{t) ~ i{tw) are equilibrated at the working 
temperature, while regions with i{t) > £{tw) are still not at equilibrium. 

Interestingly enough, we demonstrated that ordering or disordering non-equilibrium 
dynamics following a quench from higher temperature or a heating process from a 
lower temperature are characterized by a higher or lower effective temperature than 
the one of the bath. This result is consistent with the intuitive interpretation of the 
effective temperature with higher (lower) values associated to more (less) disordered 
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configurations than the equilibrium ones at the working temperature. A similar 
dependence on the initial condition was derived by Berthier et al in the 2d XY model |25j . 

The two-times length scale / also shows the latter property. In the aging regime, 
for fixed the length grows with At for all T ^ Tq while for fixed At it grows with t^, 
when To > T and it decreases with ty^ when Tq < T. The former behaviour is similar 
to what is observed in conventional glassy systems such as the 3d Edwards-Anderson 
spin glass |19] and models of particles in interaction [18j. The heating procedure was 
not studied in such cases. 

The wave-vector dependent correlation (Cg) that plays the role of the incoherent 
scattering function in studies of glassy systems is particularly interesting. We showed 
that, although {Cg) is simply related to the roughness, the characteristic multiphcative 
scaling is not easily detected in {Cg). This might be the case in other systems, such as 
colloidal glasses where it was recently shown that diffusive correlations clearly display 
multiphcative aging scaling [31] , while this was not previously reckoned in the study of 
the incoherent scattering function [T7] . 

One can also observe that the aging behaviour of (Cg) resembles strongly the 
experimental results in laponite [IT]. In particular it was found in this system that 
the incoherent scattering function displays a waiting-time dependent plateau at long 
time-delay. This suggests that a similar equilibration mechanism might be at work in 
the relaxation dynamics of the rather complex laponite samples, where some competing 
length scale is confining the particle fluctuations, thus leading to saturation of the 
incoherent scatering function. The waiting-time dependent saturation is in line with the 
fact that effective temperatures, measured through the FDT, should become the bath 
temperature at fixed time-delay and sufficiently long waiting times where saturation is 
found. The measurements in [iSO] are such that the effective temperature does indeed 
tend to the bath temperature at fixed working frequency - equivalently time-delay - 
when tw is large enough, although this result remains controversial [35]. 

With this analytic study we showed that the qualitative aging dynamics of the 
vortex glass [TOl [TT] as well as the elastic line in a quenched disordered environment [HI E] 
is mainly due to the non-equilibrium relaxation of the pure elastic lines. The effect of 
quenched disorder and line-line interactions is to change the details of the scaling, more 
precisely the temperature and time- dependence in i(t), but not the qualitative features. 
Along this line, the study of the aging dynamics of the pure KPZ equation [36] will allow 
one to better rationalize the results in [H], where the aging dynamics of this equation 
with a disordered potential - and driving force - was analyzed. Furthermore, the results 
obtained here could be strongly related to the non-equilibrium relaxation dynamics of 
confined polymers [37] . 

We presented the first analytic calculation of finite-size fluctuations during an out- 
of-equilibrium relaxation. This study complements the analysis in [33] and [23] for 
the width fluctuations at saturation and growth and in [31], [32] for other quantities 
such as the maximum height displacement - indeed, it is also simple to include the 
t^y-dependence in this calculation. Our results make explicit the dependence on the 
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waiting-time and display the crossover to equilibrium. They constitute a benchmark for 
Racz proposal to classify interface dynamics into universality classes [20] , now extending 
it to the non-equilibrium relaxation. 

As regards the time-reparametrization invariance scenario for glassy dynamics we 
do not expect it to hold, without modification, in models with multiplicative aging 
scaling. Following the steps sketched in [22] to study the asymptotic averaged dynamics 
of the EW line (or massless scalar field) and its corresponding dynamic action, one 

1/2 

soon realizes that the multiplicative t^ factor has to be extracted from the asymptotic 
analysis to search for time-reparametrization invariance. This is similar to what was 
shown in [27] for the O(A^) model in the large N limit. One should also notice that the 
dynamics of the EW line depends on the dimensionality of the transverse space. For 
instance, one can show that for infinite systems diffusion disappears and the aging 
regime persists at infinite waiting-times {i.e. the scaling becomes additive) in two 
transverse dimensions [23]. One has then the interesting possibility of testing the time- 
reparametrization invariance scenario in the 'critical' EW equation with two transverse 
dimensions. The detailed analysis of the djTianiic symmetries of the generic EW line 
goes beyond the scope of this article. 

We conclude with a note on the relevance of our results for coarsening phenomena. 
The domain walls between equilibrated regions during domain growth are usually 
described as elastic objects. Recently, the distribution of domain sizes and perimeter 
lengths in two-dimensional Ising ferromagnetic growth was shown to be unexpectedly 
non-trivial [3H] • The wide distribution of domain sizes and domain wall lengths combined 
with the highly non-trivial fluctuating dynamics of finite-length elastic lines derived 
here suggest that characterising the fluctuations of standard two-times observables in 
domain-growth phenomena can be a quite complicated problem. 

Acknowledgment s 

We thank the Universidad Nacional de Mar del Plata, Argentina, for hospitality during 
the preparation of this work and C. Chamon, D. Dommguez, T. Giamarchi, G. Schehr 
and H. Yoshino for very useful discussions. LFC acknowledges financial support from 
Secyt-ECOS R AOIEOI and PICS 3172, SB from the Swiss National Science Foundation 
under MaNEP and Division II, and JLI from CONCIET PIP05-5648 and ANPCYT 
PICT04-20075. LFC is a member of Institut Universitaire de France. 



[1] Barabasi A-L and Stanley HE, 1995 Fractal concepts in surface growth (Cambridge: Cambridge 
University Press) 
Halpin-Healey T and Zhang Y-C, 1995 Phys. Rep. 254 215 
[2] Bray AJ, 1994 Adv. Phys. 43 357 



Out- of- equilibrium Edwards- Wilkinson equation 



32 



[3] Blatter G, Feigel'man MV, Geshkenbein VB, Larkin AI and Vinokur VM, 1994 Rev. Mod. Phys. 
66 1125 

Nattermann T and Scheidl S, 2000 Adv. Phys. 49 607 
[4] Hansen A, Hinrichsen EL and Roux S, 1991 Phys. Rev. Lett. 66 2476 

Bouchaud E, 1997 J. Phys.: Condens. Matter 9 4319 

Alava M, Nukalaz PKVV and Zapperi S, 2006 Adv. Phys. 55 349 
[5] Sahimi M, 1995 Flow and Transport in Porous Media and Fractured Rock (New York: John Wiley 
& Sons) 

Alava M, Dube M and Rost M, 2004 Adv. Phys. 53 83 
[6] Edwards SF and Wilkinson DR, 1982 Proc. R. Soc. London, Ser. A 381 17 

[7] Cugliandolo LF, 2004 Slow Relaxations and Nonequilihrium Dynamics in Condensed Matter (Les 
Houches-Ecole d'Ete de Physique Theorique vol 77), ed J-L Barrat et al. (Berlin: Springer) Also 
available as |cond-mat/0210312 
[8] Yoshino H, 1996 J. Phys. A: Math. Gen. 29 1421 
Yoshino H, 1998 Phys. Rev. Lett. 81 1493 
Barrat A, 1997 Phys. Rev. E 55 5651 
[9] Bustingorry S, Iguain JL, Chamon S, Cugliandolo LF and Domi'nguez D, 2006 Europhys. Lett. 76 
856 

[10] Bustingorry S, Cugliandolo LF and Dommguez D, 2006 Phys. Rev. Lett. 96 027001 
Bustingorry S, Cugliandolo LF and Dommguez D, 2007 Phys. Rev. B 75 024506 

[11] Schchr G and Rieger H, 2005 Phys. Rev. B 71 184202 
Schehr G and Le Doussal P, 2005 Europhys. Lett. 71 290 

[12] Portier F, Kriza G, Sas B, Kiss LF, Pethes I, Vad K, Keszei B and Williams FIB, 2002 Phys. Rev. 
B 66 140511 

Exartier R and Cugliandolo LF, 2002 Phys. Rev. B 66, 012517 

Du X, Li G, Andrei EY, Greenblatt M and Shuk P, 2007 Nature Phys. 3, 111 
[13] Kolton A, Rosso A and Giamarchi T, 2005 Phys. Rev. Lett. 95 180604 
[14] Ramasco JJ, Lopez JM and Rodriguez MA, 2006 Europhys. Lett. 76 554 
[15] Kardar M, Parisi G and Zhang YC, 1986 Phys. Rev. Lett. 56 889 
[16] Cugliandolo LF and Le Doussal P, 1996 Phys. Rev. E 53 1525 

Cugliandolo LF, Kurchan J and Le Doussal P, 1996 Phys. Rev. Lett. 76 2390 

Konkoh Z, Hertz J and Franz S, 2001 Phys. Rev. E 64 051910 

Konkoh Z and Hertz J, 2003 Phys. Rev. E 67 051915 

Goldschmidt YY, 2006 Phys. Rev. E 74 021804 
[17] Bonn D, Tanaka H, Wegdam G, KeUay H and Meunier J, 1999 Europhys. Lett. 45 52 

Tanaka H, Jabbari-Farouji S, Meunier J and Bonn D, 2005 Phys. Rev. E 71 021402 
[18] Parisi G, 1999 J. Phys. Chem. 103, 4128 

Parsaeian A and Castillo HE, 2007 Nature Phys. 3 26 
[19] Jaubert LDC, Chamon C, Cughandolo LF and Picco M, 2007 J. Stat. Mech. P05001 
[20] Racz Z, 2003 SPIE Proceedings 5112 248 

[21] Chamon C, Kennett MP, Castillo HE and Cughandolo LF, 2002 Phys. Rev. Lett. 89 217201 
Castillo HE, Chamon C, Cughandolo LF and Kennett MP, 2002 Phys. Rev. Lett. 88 237201 
Castillo HE, Chamon C, Cughandolo LF, Iguain JL and Kennett MP, 2003 Phys. Rev. B 68 134442 
Chamon C, Charbonneau P, Cugliandolo LF, Reichman D and Sellitto M, 2004 J. Chem. Phys. 
121 10120 

[22] Chamon C and Cugliandolo LF, 2007 Fluctuations in glassy systems Preprmt larXiv:0704.0684[ 

[23] Antal T and Racz Z, 1996 Phys. Rev. E 54 2256 

[24] Cughandolo LF, Kurchan J and Parisi G, 1994 J. Phys. I 4 1641 

[25] Berthier L, Holdsworth PWC and Sellitto M, 2001 J. Phys. A: Math. Cen. 34 1805 

[26] Ciuchi S and de Pasquale F, 1988 Nucl. Phys. B 300 31 

Cugliandolo LF and Dean DS, 1995 J. Phys. A: Math. Cen. 28 L453 



Out- of- equilibrium Edwards- Wilkinson equation 



33 



Cugliandolo LF and Dean DS, 1995 J. Phys. A: Math. Gen. 28 4213 
[27] Chamon C, Cugliandolo LF, Yoshino H, 2006 J. Stat. Mech. P01006 
[28] Yang H-N, Lu T-M and Wang G-C, 1992 Phys. Rev. Lett. 68 2612 

Yang H-N, Lu T-M and Wang G-C, 1993 Phys. Rev. B 47 3911 
[29] Family F and Vicsek T, 1985 J. Phys. A: Math. Gen. 18 L75 
[30] Abou B and Gallet F, 2004 Phys. Rev. Lett. 93 160603 
[31] Majumdar SN and Comtet A, 2004 Phys. Rev. Lett. 92 225501 

Majumdar SN and Comtet A, 2005 J. Stat. Phys. 119 777 
[32] Schehr G and Majumdar SN, 2006 Phys. Rev. E 73 056103 

[33] Foltin G, Oerding K, Racz Z, Workman RL and Zia RKP, 1994 Phys. Rev. E 50 R639 

Plischke M, Racz Z and Zia RKP, 1994 Phys. Rev. E 50 3589 

Racz Z and Plischke M, 1994 Phys. Rev. E 50 3530 

Bramwell ST, Holdsworth PCW and Pinton JF, 1998 Nature 396 552 
[34] Wang P, Song C and Makse HA, 2006 Nature Phys. 2 526 
[35] Greinert N, Wood T and Bartlctt P, 2006 Phys. Rev. Lett. 97 265702 

Jabbari-Farouji S, Mizuno D, Atakhorrami M, MacKintosh EC, Schmidt CF, Eiser E, Wegdam 
GH and Bonn D, 2007 Phys. Rev. Lett. 98 108302 
[36] Bustingorry S, in preparation 

[37] Rahmani A, Castelnovo C, Schmit J and Chamon C Prevrint larXiv:0704. 16631 
[38] Arenzon JJ, Bray AJ, CugUandolo LF and SiciUa A, 2007 Phys. Rev. Lett. 98 145701 
Sicilia A, Arenzon JJ, Bray AJ, Cugliandolo LF, in preparation 



